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Abstract. Propagation of large-amplitude waves in plasmas 
is subject to several sources of nonlinearity due to relativistic 
effects, either when particle quiver velocities in the wave 
held are large, or when thermal velocities are large due to 
relativistic temperatures. Wave propagation in these condi¬ 
tions has been studied for decades, due to its interest in sev¬ 
eral contexts such as pulsar emission models, laser-plasma 
interaction, and extragalactic jets. 

For large-amplitude circularly polarized waves propagat¬ 
ing along a constant magnetic held, an exact solution of the 
huid equations can be found for relativistic temperatures. Re¬ 
lativistic thermal effects produce: (a) a decrease in the ef¬ 
fective plasma frequency (thus, waves in the electromagnetic 
branch can propagate for lower frequencies than in the cold 
case); and (b) a decrease in the upper frequency cutoff for 
the Alfven branch (thus, Alfven waves are conhned to a fre¬ 
quency range that is narrower than in the cold case). It is also 
found that the Alfven speed decreases with temperature, be¬ 
ing zero for inhnite temperature. 

We have also studied the same system, but based on the 
relativistic Vlasov equation, to include thermal effects along 
the direction of propagation. It turns out that kinetic and huid 
results are qualitatively consistent, with several quantitative 
differences. Regarding the electromagnetic branch, the effec¬ 
tive plasma frequency is always larger in the kinetic model. 
Thus, kinetic effects reduce the transparency of the plasma. 
As to the Alfven branch, there is a critical, nonzero value of 


the temperature at which the Alfven speed is zero. For tem¬ 
peratures above this critical value, the Alfven branch is sup¬ 
pressed; however, if the background magnetic held increases, 
then Alfven waves can propagate for larger temperatures. 

There are at least two ways in which the above results can 
be improved. First, nonlinear decays of the electromagnetic 
wave have been neglected; second, the kinetic treatment con¬ 
siders thermal effects only along the direction of propagation. 
We have approached the hrst subject by studying the para¬ 
metric decays of the exact wave solution found in the context 
of huid theory. The dispersion relation of the decays has been 
solved, showing several resonant and nonresonant instabili¬ 
ties whose dependence on the wave amplitude and plasma 
temperature has been studied systematically. Regarding the 
second subject, we are currently performing numerical 1-D 
particle in cell simulations, a work that is still in progress, al¬ 
though preliminary results are consistent with the analytical 
ones. 


1 Introduction 

Due to their interest both in astrophysical and laboratory 
plasmas, relativistic electron-positron plasmas have been a 

subject of study for decades (Rufhni et al., 2010). In astro¬ 
physics, they are relevant in systems such as accretion disks 

around black holes (Bjornsson et al., 1996; Liang, 1979; 
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White and Lightman, 1989), models of the early universe 
(Gibbons et al., 1985; Tajima and Taniuti, 1990; Tatsuno 
et al., 2003; Lesch and Pohl, 1992), supernova remnants and 
active galactic nuclei (Hardy and Thoma, 2000; Reynolds 
et al., 1996), pulsar magnetospheres (Curtis, 1991; Istomin 
and Sobyanin, 2007; Manchester and Taylor, 1977; Stur- 
rock, 1971), magnetars (neutron stars with magnetic fields 
up to ~ 10 14 G) (Beskin et al., 1993), hypothetical quark 
stars (Usov, 1998), and gamma-ray bursts (Piran, 1999, 
2004). Regarding laboratory plasmas, they have been con¬ 
sidered in the study of ultra-intense lasers (Blaschke et al., 
2006), and in laboratory and tokamak plasmas (Zank and 
Greaves, 1995). For instance, recent experiments on rela¬ 
tivistic electron-positron creation with short ultra-intense 
laser pulses (~ 10 20 Wcm -2 ) have been performed (Chen 
et al., 2009), where measurements indicate the positron den¬ 
sity to be ~ 10 16 cm -3 . Electron-positron pair production 
is also expected to occur in large tokamaks, where up to 
about ~ 10 14 positrons may be created in collisions between 
runaway electrons and background particles (Helander and 
Ward, 2003). 

Of particular interest are effects related to wave propaga¬ 
tion in these plasmas, such as proposed pulsar radio emission 
processes (Luo et al., 2002), bulk acceleration of relativistic 
jets (Iwamoto and Takahara, 2002), jet formation (Sawyer, 
2008; Wardle et al., 1998), electron-positron pair annihila¬ 
tion into one photon in the presence of a strong magnetic 
field (Harding, 1986), and also in laboratory environments, in 
problems such as pair production by optical lasers (Blaschke 
et al., 2006; Chen et al., 2011). 

In several of the environments mentioned above, relati¬ 
vistic effects and temperature play an important role, thus 
it is fundamental to understand wave propagation modes in 
relativistic plasmas with temperature. In this paper we will 
focus on the particular case of circularly polarized electro¬ 
magnetic waves, which, though simple, allows us to study in 
detail the effect of relativistic temperatures on wave propa¬ 
gation in relativistic hot plasmas. 

First, in Sect. 2 it is shown that an exact solution can 
be found for the fully relativistic fluid equations, including 
relativistic temperatures (Asenjo et al., 2009). The disper¬ 
sion properties of the wave are studied, showing that two 
branches appear, an Alfven branch and an electromagnetic 
one, and also showing that relativistic effects increase the 
plasma transparency, and confine the Alfven branch to a 
lower frequency and wave number range. Then, in Sect. 3 
the problem is considered from the point of view of a ki¬ 
netic model (Dominguez et al., 2012). Although qualitative 
agreement with the fluid results is found, a major difference 
between both models is the fact that, in the kinetic one, the 
Alfven branch is suppressed if the plasma temperature is 
large enough. 

There are at least two ways in which the above results can 
be improved. First, nonlinear decays of the electromagnetic 
wave have been neglected; second, the kinetic treatment con¬ 


siders thermal effects only along the direction of propagation. 
In order to address the first subject, we have also studied the 
parametric decays of the exact wave solution found in the 
context of fluid theory (Lopez et al., 2012), results that are 
shown in Sect. 4. There are several resonant and nonreso¬ 
nant instabilities, whose dependence on the wave amplitude 
and plasma temperature are studied systematically. Regard¬ 
ing the second subject, thermal effects in all directions with 
respect to the background magnetic field can be readily taken 
into account in numerical simulations. Thus, we are currently 
performing numerical 1-D particle in cell simulations, a work 
that is still in progress, (Lopez et al., 2014) although prelim¬ 
inary results, shown in Sect. 5, are consistent with the ana¬ 
lytical ones. Finally, in Sect. 6 results are summarized and 
discussed. 

2 Exact solution of relativistic fluid equations 

Exact solutions for the plasma equations can be found for 
cold nonrelativistic plasmas. For instance, circularly pola¬ 
rized Alfven electromagnetic waves propagating parallel to 
an external magnetic field are an exact solution of the mag¬ 
netohydrodynamic equations even when the amplitude is 
large (Barnes and Hollweg, 1974). Also, a circularly po¬ 
larized wave in a multiple ion species plasma with drifts 
is a finite-amplitude solution of the cold plasma model 
(Gomberoff et al., 1994). The nonlinear propagation of cir¬ 
cularly polarized electromagnetic waves in unmagnetized 
electron-positron-ion plasmas has also been studied in the 
cold (Berezhiani and Mahajan, 1994) and relativistically hot 
(Berezhiani and Mahajan, 1995) cases, showing the existence 
of stable localized structures. 

Here we propose an approach that permits one to find 
an exact solution for the propagation modes in a relativistic 
electron-positron plasma with constant, arbitrary tempera¬ 
ture, within the context of a fluid theory. This can be done 
by basing our approach on the magnetofluid field unification 
formalism (Mahajan, 2003). 

In this unification approach, the whole plasma is treated as 
a unique field where the electromagnetic field is coupled with 
the charged fluid field through a function that carries statisti¬ 
cal information of the system. This leads to a simple and el¬ 
egant way of describing relativistic plasmas. The formalism 
has also been used successfully to study wave propagation in 
relativistic two-fluid plasmas (Soto-Chavez et al., 2010), to 
find equilibrium states via a variational principle (Pino et al., 
2010), and has been extended to non-Abelian fields in order 
to study quark-gluon plasmas (Bambah et al., 2006). 

In this paper, we use the approach presented by Mahajan 
(2003) to derive the dispersion relation for circularly pola¬ 
rized waves of arbitrary amplitude propagating along a con¬ 
stant magnetic field, for arbitrary temperatures. 
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2.1 Magnetofluid unification 


the continuity equation 


Usually, the interaction of particles with electromagnetic 
fields is described by introducing a “minimal coupling” in 
the momentum, p —> p — e A/ c, where A is the vector poten¬ 
tial. This leads to an energy-momentum conservation equa¬ 
tion. This equation and Maxwell equations describe the basic 
dynamics of charged relativistic particles in plasmas. 

It has been suggested (Mahajan, 2003) that the coupling 
of a relativistic charged fluid, at a given temperature, with 
the electromagnetic field can be described by an antisym¬ 
metric field tensor that contains the statistical information of 
the system. Thus, the set of equations to describe the plasma 
dynamics in the homentropic regime, for species /, are the 
Maxwell equations, the continuity equation, and the equation 

<ljUj»M v /=0. ( 1 ) 

The field M^ v is the tensor that couples the electromag¬ 
netic and fluid fields. This tensor is defined as M^ v = F^ v + 
( mjC 2 /qj ) S ^ v , with F^ v as the electromagnetic tensor, qj is 
the charge, and where the new antisymmetric tensor S f j V = 

3 M — 3 V is introduced, representing the re¬ 

lativistic thermal fluid. Here, Uj' —> (yj, YjVj/c ) is the four- 
velocity of the fluid, vj is the species velocity, 



and c is the speed of light. 

The parameter fj is a function of the temperature T. An 
explicit form for f(Tj) can be obtained by assuming a given 
statistical behavior for the gas. For instance, if the system fol¬ 
lows a Maxwell-Juttner equilibrium distribution (Berezhiani 
and Mahajan, 1995; Mahajan, 2003), 



where K 2 and K 3 are the modified Bessel functions of or¬ 
ders 2 and 3, respectively, and ks is the Boltzmann constant. 
However, within the treatment carried out in this paper, no 
explicit description for the function / is needed, and all the 
analytical and numerical results that follow are independent 
of the underlying particle distribution function. 

Then the plasma dynamics turns out to be given by the 
momentum equation (Mahajan, 2003) 

(s + -'- v )= |( E V' x, )-i v "' <4) 

which is the spacelike component of Eq. (1), with pj the 
pressure of species j, the pressure equation 

(It + v ‘ ■ v ) Pi = mi n c * (h + ■ v ) fi • (5) 


i^ + V.( njVj) = 0 , ( 6 ) 

and Maxwell equations. These yield 
V • E = 4j re (n v — n e ) > (7) 


4 7 T 13 E 

B = - J + -—, 
c c at 


and, in the Lorentz gauge, 


1 3 2 A 

c 2 dt 2 


-V 2 A = 


4tt 471 

— J = — 

c c 


Yww • 

j 


( 8 ) 


( 9 ) 


Also, we will consider pj = HjTj, where nj = nj /yj is 
the plasma density in the rest frame (Asenjo et al., 2009). 

An alternative form of Eq. (4) can also be found in terms 
of the enthalpy density h (Gomberoff and Galvao, 1997). 

The formalism of magnetofluid unification (Mahajan, 
2003) is very general and allows us to study the dynamics of 
a charged fluid (plasma) in an electromagnetic field in a uni¬ 
fied way, treating them as a single field. On the other hand, 
it also provides a general framework to study various plasma 
physics phenomena, taking into account relativistic temper¬ 
ature effects in a consistent way. In particular, the effect of 
relativistic temperatures on wave propagation in plasmas can 
be studied systematically. For the sake of simplicity, here we 
will consider an equal mass electron-positron plasma, where 
a circularly polarized electromagnetic wave propagates along 
a constant background magnetic field. 


2.2 Circularly polarized waves 


We consider a circularly polarized wave that propagates in 
the direction of a constant background magnetic field along 
the z axis. Thus, the electric and magnetic fields are given by 

Eo(z, t) = E [sin (koz — (Dot)x -cos(koz - <z>o t)y] , (10) 

Bo(z,t ) = B [cos (koz - mt) x + sin (koz - coot) y] + Bq z z, (11) 


respectively, where ko is the wave number and coo is the fre¬ 
quency of the monochromatic wave. 

The amplitudes of the electric and magnetic fields can be 
related by Maxwell’s equations as 



( 12 ) 


Electrons and positrons have equal constant densities n v = 
n e = n, and equal constant drift velocities v z o = v z oz, so that 
their velocity can be written as Vj = v t j+v z 0 , where v t j is 
transverse with respect to the background magnetic field. We 
will also denote their mass by m p = ra e = m, and their charge 
by <7p = = e. 
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In principle, we assume that electrons and positrons have 
constant but different temperatures, so that / e # / p . It turns 
out, as we will see below, that the purely transverse circu¬ 
larly polarized wave is an exact solution of the field equa¬ 
tions, which is consistent with the assumption of no pressure 
fluctuations, that is, no fluctuations in /. 

The particle velocities induced by the wave field are purely 
transverse, hence the amplitude of the circularly polarized 
velocity is proportional to the constant amplitude of the cir¬ 
cularly polarized vector potential field (which can be easily 
shown in the Lorentz gauge, for instance), and the relativistic 
factor, which can be written as 



is constant for both species. 

This suggests that it is convenient to rewrite Maxwell’s 
equations in terms of the transverse velocity v t j, yielding a 
wave equation for the vector potential 



1 9 2 \ 

C 2 dt 2 ) 


A 



(14) 


where 


j'-T.wj ( 15 > 

j 

is the transverse current. The wave Eq. (14), Poisson’s equa¬ 
tion for the electrostatic potential 

W 2 cp = —47r Ygnj , (16) 

j 


the motion Eq. (4), and the continuity equation 
dn j / n. 

-^- = -V.(njv j ), (17) 

are the complete set of fluid equations that describe the 
plasma. 

Defining complex transverse quantities in the form 


C± = C x +iCy, (18) 

we find that 

Cl = Ce> (koz -aot ). (19) 


Then, the particle velocity in the wave field can be written in 
terms of the wave magnetic field as (Asenjo et al., 2009) 


co 


qjB ± 


fjYj(o f — Q C j mck 


( 20 ) 


where Q c j =qj #o /me and co' = co — kv z o are the gyrofre- 
quency for both species and the Doppler shifted frequency, 
respectively. 


Finally, this leads to the following dispersion relation for 
an electron-positron plasma (Asenjo et al., 2009): 



( 21 ) 


with co v = y/47tne 2 /m the electron plasma frequency. This 
dispersion relation is an exact solution of the plasma equa¬ 
tions for finite-amplitude circularly polarized propagating 
waves along the magnetic field in a relativistic plasma with 
temperature. 

It is interesting to note that Eq. (21) differs from the cold 
plasma case (Matsukiyo and Hada, 2003) only by a fac¬ 
tor fj. We can understand this as follows. For single par¬ 
ticle motion, it is possible to convert a non-relativistic re¬ 
sult into a relativistic one by changing rrij —> Yj m j- This is 
possible because mass is involved in the momentum equa¬ 
tion. However, this simple replacement is not possible for 
a fluid, since velocity is related nonlinearly to momentum 
through Yj- Therefore, the average on momentum (pj) is not 
equivalent to (Yj)(mjV j). We can think, though, that there 
is some proportionality factor between both quantities, say 
fj, such that (pj) = fj{Yj){rrijVj). Thus, it would be pos¬ 
sible to convert a nonrelativistic fluid force equation into a 
relativistic one with the same prescription as for single par¬ 
ticles, ( ntjVj) -> (Pj), but, in order to take into account the 
statistics, rrij — > fjYj m j • This is exactly what is needed to 
go from the cold fluid dispersion relation (fj = Yj = 1 in 
Eq. 21) to Eq. (21). 

For / p = / e = 1 and no drift, the dispersion relation for a 
cold relativistic plasma is recovered (Matsukiyo and Hada, 
2003; Gomberoff and Galvao, 1997). The same result as in 
Matsukiyo and Hada (2003) can be obtained from a kinetic 
approach (Munoz et al., 2006), in the limit where the velo¬ 
city distribution is a Dirac delta (see, e.g., Eq. (23) in Munoz 
et al., 2006). 

At the opposite limit of very high temperature, the trans¬ 
verse mode becomes a light mode c 2 k 2 = co 2 , consistent with 
the relativistic decrease in the effective plasma frequency. 


2.3 Analysis of the dispersion relation 


In the rest of this paper, we will consider the case where both 
species have the same temperature, so that / p = / e = /. In 
order to study the dispersion relation Eq. (21) we will nor¬ 
malize all frequencies to the positron gyrofrequency Q cp = 
Q c , and velocities to the speed of light c. It is convenient 
to define two adimensional parameters, related to the plasma 
frequency and to the wave amplitude, 


g|A| _ e\E ± \ _ g|fl|_| 
Q 2 ’ me 2 me co mc 2 k 


( 22 ) 


Notice that a corresponds to the particle transverse momen¬ 
tum due to the wave. 
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In order to plot the dispersion relation, Eqs. (13), (20), 
and (21) are solved simultaneously for yj, Vj±, and co , for 
a given k , as outlined in Munoz et al. (2006). 

In Fig. 1 the dispersion relation (21) is plotted for various 
values of /, and for a = 1, a = 0.1. There are two branches: 
an electromagnetic branch, with a lower cutoff at the effec¬ 
tive plasma frequency, and an Alfven branch, which has an 
upper cutoff in wave number, and an upper cutoff in fre¬ 
quency. 

The solid line corresponds to the relativistic cold plasma, 
/ = 1. It can be seen that, as mentioned in Sect. 2.2, when 
the temperature increases the effective plasma frequency de¬ 
creases, and the electromagnetic wave becomes a light wave. 
This effect can be better appreciated in Fig. 2, where the ef¬ 
fective plasma frequency <z> p ff (given by the lower frequency 
cutoff for the electromagnetic branch in Fig. 1) is plotted for 
the same values of / used in Fig. 1. <z> p ff also decreases due 
to the relativistic effect on the mass, and thus we plot it as 
a function of a as well. However, for high enough tempera¬ 
tures, the variation of the plasma frequency with wave am¬ 
plitude is negligible. 

Regarding the Alfven branch, Fig. 1 shows that it starts 
at the origin following the usual linear dispersion relation, 
co = VAk. In this region, y p , y e ~ 1. Then, the Alfven velo¬ 
city can be obtained from Eq. (21) by first rewriting it in the 
form: 


c 2 k 2 = co 2 — 2 fco 2 


co 


(fco) 2 


■al 


If cof « £2 C , which is always satisfied for low enough 
frequencies, we find that the Alfven velocity is given by 


v A = 


c 

J1+2 f 


(23) 


As wave number increases along the Alfven branch, there 
is an upper frequency cutoff <z> cr it given by 


fr^crit 

£2 C 


1 

7 




(24) 


Thus, when temperature increases, Alfven waves are eventu¬ 
ally confined to a very narrow frequency band. This can also 
be seen in Fig. 1. 

Another interesting feature of the Alfven branch is the 
existence of a maximum value for the wave number, k max . 
This upper cutoff occurs because one species (positrons, in 
this case) becomes ultra-relativistic, essentially moving with 
the speed of light in the wave field, thus, y p = oo. Electrons 
do not resonate with the wave, so y e is finite. In order to 
obtain an analytic expression for the wave number cutoff for 
the Alfven branch, let us first notice that, from Eq. (20), it 



Fig. 1. Electromagnetic and Alfven branches for the general disper¬ 
sion relation, Eq. (21), for a = 1, a = 0.1. Solid line: cold plasma 
case (/ = 1); dashed line: f = 2; dotted line: / = 10. 



Fig. 2. Effective plasma frequency &> p ff as a function of wave am¬ 
plitude a, for various temperatures. Solid line: cold plasma case 
(/ = 1); dashed line: f = 2; dotted line: / = 10. 


follows that 


fyjco-Qcj a c 


(25) 


where the plus (minus) sign corresponds to positrons (elec¬ 
trons), so the dispersion relation Eq. (21) can be written as 


60 " 


C 2 k 2 = CO 2 - ~(v± p - V ± e ) • 

ac x 


(26) 


The cutoff occurs for co = 0, and in this limit v± p ~ — c (no¬ 
tice that the left-hand term in Eq. 25 is negative along the 
Alfven branch for all values of k ), whereas v± e becomes neg¬ 
ligible. Thus, it follows that the maximum wave number is 
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given by 

(O v 

fcmax = —£= ■ (27) 

c^ya 

This is consistent with the numerical result in Fig. 1. Notice, 
in particular, that & max depends only on wave amplitude, not 
on temperature. One could argue that, if the transverse par¬ 
ticle velocity is ultra-relativistic, then thermal velocities are 
not relevant. 

Considering the change of sign of dco/dk in the Alfven 
branch, for k < &(<z> C rit) and k > k(co CY it), we will identify 
these regions as the normal and the anomalous zone of the 
Alfven branch, respectively. 

3 Kinetic model 


electron-positron plasma is solved. This corresponds to a 
generalization of the numerical results shown in Sect. 2 and 
by Asenjo et al. (2009), for the case of a plasma with con¬ 
stant longitudinal drift velocity. Then, the basic equations 
of the kinetic model are presented in Sect. 3.2. In Sect. 3.3 
the dispersion relation of this wave is calculated. Then, in 
Sect. 3.4 the kinetic dispersion relation obtained in Sect. 3.3 
is solved numerically, and results are compared with the cor¬ 
responding fluid treatment. These results can also be found 
in Dominguez et al. (2012). 

3.1 Fluid dispersion relation with drift 

We start by defining normalized variables x = co/Q c , y = 
ck/Q c, and /3 = p z /(mc), so that Eq. (21) can be written in 
normalized form as: 


In Sect. 2, an exact solution in the context of fluid theory was 
obtained for a circularly polarized wave propagating along a 
constant magnetic field in an electron-positron plasma with 
temperature, and its temperature dependence was studied nu¬ 
merically. Relativistic thermal effects produce: (a) a decrease 
in the effective plasma frequency (thus, waves in the electro¬ 
magnetic branch can propagate for lower frequencies than in 
the cold case); and (b) a decrease in the upper frequency cut¬ 
off for the Alfven branch (thus, Alfven waves are confined to 
a frequency range that is narrower than in the cold case). 

These results have been found in the context of a fluid 
model. However, effects such as Landau and cyclotron damp¬ 
ing are absent in such a model, and thus it would be interest¬ 
ing to study to what extent the fluid results are modified by 
kinetic effects. Many works have been devoted to the study of 
kinetic effects on wave propagation in relativistic electron- 
positron plasmas, either in the absence of a magnetic field 
(Munoz, 2004; Liu and Liu, 2011), to analyze the counter¬ 
streaming plasmas in pulsar environments (Verdon and Mel¬ 
rose, 2011), to study the effect of nonlinearities and wave 
damping (Chaudhary et al., 2010), to model the propagation 
of Bernstein waves (Gill and Heyl, 2009), etc. 

In this section, we study the propagation of a large- 
amplitude circularly polarized electromagnetic wave in a 
magnetized plasma, analyzing the consistency of the results 
obtained in Sect. 2 with those obtained from a kinetic treat¬ 
ment (Munoz et al., 2006). In effect, the same system is 
studied, but starting from the relativistic Vlasov equation. 
The corresponding kinetic dispersion relation is numerically 
studied for various temperatures, and results are compared 
with the purely fluid treatment. In our treatment, we assume 
that kinetic effects are only important along the propagation 
direction, which makes it necessary to generalize the results 
in Sect. 2, and to study the fluid dispersion relation Eq. (21), 
but in the presence of drift. 

In order to gain intuition for developing the posterior ki¬ 
netic analysis, in Sect. 3.1, the fluid relativistic dispersion 
relation of a circularly polarized electromagnetic wave prop¬ 
agating along a constant background magnetic field in an 


0 = y 


2 



xyj - yP i 

/ (x Yj - v/3) - <Jj Yj ’ 


(28) 


where aj is the sign of the charge of species j, so that cr p = 1 , 
and cr e = — 1. 

In order to solve the fluid dispersion relation given by 
Eq. (28), the roots of the following quartic equation for yj 
must be found: 


^ ~ Yj f 2 ~ yf~~ (yPf + °j) 


+ r 


{ypf + ojf 


— (l + £ 2 )/ 2 —< 


+ Yj - [0 + P 2 ) f (yPf+ a j) + y/te 2 ] 
-r^[( l + P 2 ) (yfi + °jf + y 2 P 2 <x 2 ] ■ 


(29) 


Equation (29) is the generalization of Eq. (15) in Munoz 
et al. (2006) for / ^ 1 and ^ 0. In Sect. 2 the case f ^ 1 
was solved (Asenjo et al., 2009). We now consider the gen¬ 
eral case given by Eq. (29), to solve the dispersion relation 
given by Eq. (28) for both nonzero temperatures (/ ^ 1) and 
nonzero drift velocities (/3 ^ 0). 

To do this, we first choose fixed values for y, /3, / and 
a , and a trial value for v. Equation (29) is then solved for 
]/ p and y e . Only solutions yj > 1 are physically acceptable. 
This yields several possible pairs (y p , y e ). Each possible pair 
corresponds to a different branch of the dispersion relation. 
When a pair is chosen, then the right-hand size of Eq. (28) 
can be evaluated. However, it turns out that, in general, the 
number of acceptable solutions of Eq. (29) for ]/ e and y v is 
greater than 2, and this number changes for each value of 
This issue complicates the choice of yj . 
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This is illustrated in Fig. 3. For given values of y, 
a , and /, we solve Eq. (29) for a certain value of x 
and /3, obtaining four roots, y> . We sort them such that 

Rey^ <Reyj 2 ^ <Rey*p <Re Figure 3 shows, for a cer¬ 
tain range of values of x and ft, the regions where the sth 
branch for positrons, y p s \ is physically acceptable, that is, 
real and greater than 1. 

Figure 3 shows that for x > 0.8, approximately, there is 

(3) 

only one possible choice, y p . However, for lower values 
of the frequency, there is a range of values of ft where 
more than one choice of y p is possible. For instance, for 
x = 0.400087, y = 0.3 and /3 = 4.4, the three possible values 
are y p (1) = 4.12392, y p (2) = 5.72992, and y p (3) = 5.86722. 

With a similar analysis in the case of y e , now we choose 
one of the possible branches for y p , and one of the branches 
for y e , and we can calculate the right-hand side of Eq. (28) 
by solving for the value of x for which it is zero. Notice that 
it could still be the case that a chosen pair (y p , y e ) yields an 
Eq. (28) with no real solutions for x. It turns out that there is 
only one choice of the (y e , y p ) pair for which a real root of 
Eq. (28) can be found. 

In Fig. 4 we plot the solutions of the dispersion rela¬ 
tion (28), obtained by the method outlined above. The solid 
line, /3 = 0, corresponds to the solution previously shown in 
Fig. 1 (Asenjo et al., 2009), which shows that our treatment 
leads to the correct limit in the absence of drift. 

Notice that the Alfven branch is much more sensitive to 
variations in drift than the electromagnetic branch. For small 
values of the frequency, it follows the usual linear dispersion 
relation, co~VAk, with va the Alfven speed, va can be cal¬ 
culated by noting that, in this region, y ± p , y± e ~ 1 + /3 2 , and 
by setting co Q c in Eq. (21). This yields 


V A 2fPy/l+P + V A0 J(l+p2)v 2 A0 + 2fy/l+fP 

C ~ *4x/ 1 +^ 2 + 2/( 1 + /S 2 ) 

where 

VA0 = - 

<Z>p 


(30) 


(31) 


is the Alfven velocity in the nonrelativistic cold fluid case. 

For /3 = 0, the result in Eq. (23) is recovered, which we 
can rewrite as 
va vao 

. (32) 


; A0 


+ 2 / 





In Fig. 5 we plot the Alfven velocity. In principle, v a depends 
both on the drift velocity and the temperature. For illustration 
purposes, in Fig. 5 we have first taken a fixed temperature 
(Fig. 5a), and then a fixed drift velocity (Fig. 5b). Qualita¬ 
tively similar results are obtained for other choices of / and 
p. What is important to notice is that the drift velocity and the 
temperature have the opposite effect on the Alfven speed: va 
increases with and decreases with /. 


Fig. 3. Isocontours of real solutions of Eq. (29) as a function 
of normalized drift j3 and normalized frequency x, for positrons. 

(5) 

Each panel corresponds to a different branch, y p , for the roots 
of Eq. (29). A few contour levels are shown. Physically acceptable 
solutions are in the regions where y(x, ft) is real and greater than 
1. Fixed values y = 0.3, a = 0.1, and f — 1 were taken, (a) 5 = 1, 
(b) s = 2, and (c) s = 3. 
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y 


Fig. 4. Solutions of the dispersion relation Eq. (21) for / = 1, 
a = 0.1, vao — 1 (see Eq. (31) below), and three values of /?: = 0 

(solid line), /3 = 0.1 (dashed line), = 0.2 (dotted line). 


3.2 Kinetic model 


In our treatment, we consider kinetic effects only along a 
background magnetic field Bo z z , and a cold fluid model in 
the perpendicular direction (Munoz and Gomberoff, 2002; 
Munoz, 2004). 

As in Sect. 2, the electromagnetic fields are described by 
Eqs. (14)—(16), whereas the fluid is described by the continu¬ 
ity equation (17) and, in the transverse direction, the motion 
equation (4) with zero temperature: 

(| + „ r v)( r „,) = £L( £+ I„ jX *). (33) 

For the longitudinal direction, we consider the one¬ 
dimensional Vlasov equation: 

^i + vWgj + ^(E + -vxB)-V vg j=0, (34) 
dt Mj \ c ) 

where gj is the velocity distribution function for species j, 
normalized so that 

J du gj(r,v,t)=njo, 

with it jo as the equilibrium unperturbed density for species 
j. and Poisson Eq. (16). 

Thus, the total density and transverse current are given by 



dvgj 0r, v,t ), 



du gj(r,v,t) Vj t , 


(35) 

(36) 


respectively. 



/ 


Fig. 5. Alfven velocity, Eq. (30), for a = 0.1, v^o = 1. (a) For fixed 
temperature, f = 1, as a function of the drift velocity ft. (b) For 
fixed drift velocity, /3 = 0, as a function of temperature, parameteri¬ 
zed by /. 


3.3 Kinetic dispersion relation 

As in Sect. 2, we now consider an electromagnetic wave 
propagating in the electron-positron plasma. In the absence 
of the wave, the plasma is considered to have no electric field, 
a uniform background magnetic field Bo z z, equal densities 
for each species (n e o = n p o = no), and equal drifts for each 
species (vjo = v z oz). 

We consider a circularly polarized wave that propagates 
along the £ axis, whose electric and magnetic fields are given 
by Eqs. (10) and (11). 

In terms of the transverse quantities Eq. (18), Maxwell and 
Vlasov equations yield: 
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On the other hand, in order to introduce kinetic effects, we 
consider the following form of the distribution function: 

gj =nogj (z, Pz,t)8^p x - pl)s(p y - pf) , (38) 

where pf is the transverse particle momentum obtained with 
the cold fluid theory. In the absence of the wave, the adimen- 
sional distribution function in Eq. (38) must satisfy g e o = 
g p o, so that densities are equal. 

Then, Eq. (37) can be written in the form: 

= /d p±v j± gj(z,p z ,t). (39) 

Note that we need to know vj± in order to be able to 
solve the dispersion relation. To do this, we consider that in 
the direction perpendicular to the magnetic field, the equa¬ 
tions governing the system are the fluid equations, so that the 
transverse velocity is given by Eq. (20). In Dominguez et al. 
(2012) it is shown that the following dispersion relation is 
obtained: 

0 = k 2 c 2 - CO 2 + col E f • < 4 °) 

3.4 Numerical solution of the kinetic dispersion relation 


We now solve the kinetic dispersion relation Eq. (40). First, 
it is convenient to write it in normalized form: 


2 

c J=P,e v 


E/ 

l—rt 


dj8 


xyj - Py gjo 

XYj-py-Cj Yj 


= 0 . 


(41) 


There are two main issues to consider at this point. One 
is the choice of the distribution function, and the other is the 
choice of yj . 

Regarding the distribution function, since we are inter¬ 
ested in comparing the kinetic dispersion relation with the 
results in Sect. 2, which include the effects of relativistic tem¬ 
perature, we take g p o = g e o = go as a Maxwell-Boltzmann- 
Juttner distribution: 


go os, H) 


1 e - n ( i +£ 2 ) 1/2 . 
2mcK\(ii) 


Thus, in the kinetic treatment, thermal information is con¬ 
tained in the normalized inverse temperature parameter fi. In 
the fluid treatment described in Sect. 2, thermal information 
is contained in the parameter /. 

Once the distribution function is chosen, the numerical 
strategy is analogous to the one outlined in Sect. 3.1, ex¬ 
cept that an integration along the /3 axis must be done. This 
means that, for each trial frequency, the corresponding val¬ 
ues of y (13) must be found. Since, in the transverse direction, 
we consider a cold fluid model, the problem of finding y is 
equivalent to the one discussed in Sect. 3.1 for a cold fluid 


with a given longitudinal drift /3. Thus, the equation to solve 
in order to find y is Eq. (29), with / = 1. It must be solved 
for each point in the /3 axis, and once the integration is com¬ 
plete, the right-hand size of Eq. (40) can be evaluated. Thus, 
rather than finding a value of, say, y e to generate the electro¬ 
magnetic branch in the fluid dispersion relation, the proper 
infinite set of values y & (p) must be found. 

If only two physically acceptable values of yj were ob¬ 
tained for each value of f3, then the problem would still be 
easy, since one of them would correspond to the electro¬ 
magnetic branch, and the other to the Alfven branch, and 
there would be no ambiguity in the calculation of the inte¬ 
gral. However, as shown in Sect. 3.1, several values of yj 
are possible, both for electrons and positrons. Moreover, the 
problem is now more complicated, since in the kinetic case, 
there is an integral over j3, thus a choice of branches for yj 
must be made at each point of the integration path, and since 
the number of valid values of yj changes with (3 (see Fig. 3), 
there is no guarantee that a choice for a given is valid along 
the complete integration path. In fact, we have found that this 
is not the case in general. By studying several families of pa¬ 
rameters, we concluded that it is very difficult to find a gen¬ 
eral, reliable algorithm to automate the choice, and some trial 
and error is still involved in the results presented here. 

The rest of the numerical procedure is equivalent to the 
fluid case: once a proper choice of yj(x,/3) is made, it is 
used to find a value of x such that Eq. (41) is satisfied, and 
the process is repeated for a different value of x . 

Figure 6 shows the solution of the dispersion relation 
Eq. (21) for three different temperatures. This figure should 
be compared with Fig. 1, where three temperatures corre¬ 
sponding to / = 1, 2 and 10 were considered. 

We should first notice that, in the Alfven branch, which 
starts at the origin, there is a small-frequency region where 
no modes are found (in this case y ~ 2). This is a result 
of the technical issues regarding the choice of yj at each 
point of the integration path, and illustrated in Fig. 3. On the 
other hand, the electromagnetic branch (for all wave num¬ 
bers) and the Alfven branch near the origin (for small values 
of x = co/Q c and y = ck/ Q c ) can always be obtained with¬ 
out ambiguities in the choice of yj . 

Thus, we will focus the remainder of the discussion on the 
electromagnetic branch and the region of the Alfven branch 
near the origin. Solving the numerical issues near the maxi¬ 
mum frequency in the Alfven branch may involve either im¬ 
proving the numerical algorithms, or improving on the cold 
model for the transverse motion (Eq. 38). In fact, prelimi¬ 
nary results from PIC simulations (Sect. 5) suggest that the 
Alfven branch is damped in the anomalous zone, and we plan 
to study this elsewhere. 

The first result to notice is that, indeed, the same branches 
as in the fluid case are found, an electromagnetic and an 
Alfven branch. For very low temperatures, the cold result is 
recovered (see solid line in Fig. 4). When the temperature is 
increased (/1 decreases), it is observed that the effect of the 
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Fig. 6. Solutions of the dispersion relation Eq. (40) for a = 0.1 and 
i^O = 1. Solid line/x = 700 (/ ~ 1); dashed line /x = 50 (f — 1.1); 
dotted line /x = 3 (/ ~ 2); dot dashed line /x = 0.4 (/ ~ 10). At the 
scale of the figure, the dashed and solid lines for the electromagnetic 
branch cannot be distinguished. Also, the dotted and dot dashed 
lines for the Alfven branch cannot be seen; this occurs because, for 
high temperatures, the Alfven branch disappears, as discussed later 
in Fig. 8. 


temperature is to decrease the effective plasma frequency, 
and to decrease the upper frequency cutoff for the Alfven 
branch. It is interesting to observe that this is qualitatively 
equivalent to what is obtained in the fluid theory. 

However, several quantitative differences are found. First, 
in the kinetic model, the effective plasma frequency is less 
sensitive to variations in temperature than in the fluid model. 
This is shown in Fig. 7, where (solution of the electro¬ 
magnetic branch, Eq. 41 with y = 0) is plotted as a function 
of inverse normalized temperature /x, for both models. For 
any given temperature, the plasma frequency is always higher 
in the kinetic case than in the fluid case. 

Regarding the Alfven branch, analytical results can be ob¬ 
tained from Eq. (40), which is consistent with the numeri¬ 
cal procedures that yield Fig. 6. Taking co <^Q C in Eq. (40) 
yields the Alfven velocity 


va 

c 


v\o-h(n) 
+ h ( m ) 


where 


h(M) = 2 / dfi — , 

j yr 

h(n) = 2 J dpg 0 (P,n)y/l + fP. 


(42) 


(43) 

(44) 


The kinetic result Eq. (42), as well as the fluid result, 
Eq. (23), are plotted as a function of /x in Fig. 8, for various 
selected values of v A o- 


Fig. 7. Effective plasma frequency as a function of inverse normal¬ 
ized temperature /x. Solid line: fluid model. Dashed line: kinetic 
model. 




Fig. 8. Alfven speed as a function of inverse normalized tempera¬ 
ture of /x. (a) Eq. (42) and (b) Eq. (30) (/3 = 0), for various values 
of v A0 . 
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We notice that increasing the temperature has the effect of 
decreasing the Alfven speed (Fig. 8a). The same qualitative 
behavior occurs in the fluid case (Fig. 8b). However, in the 
kinetic treatment, unlike the fluid case, the Alfven branch is 
suppressed for high temperatures. The threshold temperature 
(proportional to /x _1 ), below which the Alfven branch ex¬ 
ists, is larger if vao increases. That is, if, say, the density is 
constant, the Alfven branch exists for a larger range of tem¬ 
peratures when the background magnetic field is increased. 
The threshold can be obtained from Eq. (42), by noting that 
it satisfies 

Vao — ^1 (Mmin) = 0 • (45) 

Figure 9 shows the value of /x m i n for various values of vao- If 
the density is fixed, Fig. 9 shows that Alfven waves can prop¬ 
agate for larger values of the temperature (/x m i n is smaller) if 
the background magnetic field is increased (vao increases). 

Figure 8 also shows that, for a given value of vao (fixed 
amplitude of the background magnetic field, for instance), 
the maximum Alfven speed is obtained in the cold limit. 

The result obtained in Eq. (42) is fully consistent with the 
fluid model in the cold limit. In fact, if /x —> oo in Eqs. (43) 
and (44), then I\ -> 0 and h —> 2, thus recovering Eq. (23). 

In Fig. 10 we plot the maximum value of the Alfven speed 
as a function of vao as given by the kinetic model (Eq. 42 
with a very large value of /x), and by the fluid model (Eq. 23, 
with / = 1). 

The fact that both results match shows that, in spite of the 
technical issues discussed above regarding the proper choice 
of y at each point in the integration path, the kinetic model 
presented here still yields the fluid result in the proper limit. 


4 Parametric decays 




Fig. 10. Maximum Alfven speed as a function of vao> obtained from 
solving the kinetic dispersion relation Eq. (40) for very low temper¬ 
ature (solid line), and from the fluid model, Eq. (23) (dotted). 


4.1 Parametric decays 

In this section we will consider the parametric perturba¬ 
tions of finite-amplitude circularly polarized electromagnetic 
waves in a relativistic electron-positron thermal plasma, 
which were found in Sect. 2. Although simple, this anal¬ 
ysis will allow us to study in detail the effect of rela¬ 
tivistic temperatures on wave propagation and its decay 
in relativistic hot plasmas. These decays have been stud¬ 
ied for electron-positron plasmas in the weakly relativistic 
limit for linear polarization (Gomberoff et al., 1997; Munoz 
and Gomberoff, 2000), and in the magnetized (Munoz and 
Gomberoff, 1998b) and nonmagnetized cases (Munoz and 
Gomberoff, 1998a) with circular polarization. Similarly, de¬ 
cays of circularly polarized Alfven waves have been studied 
(Matsukiyo and Hada, 2003), also in the weakly relativistic 
regime. On the other hand, the effect of relativistic tempera¬ 
tures has not been considered in these studies. 


In Sect. 4.3, the dispersion relation for the parametric de¬ 
cays is presented. In Sect. 4.4, this dispersion relation is 
solved numerically and several features are discussed. 

4.2 Exact solution 


In the parametric decay analysis, the exact circularly polari¬ 
zed electromagnetic wave found in Sect. 2 corresponds to the 
zeroth-order solution. Thus, zeroth-order fields are given by 
Eqs. (10) and (11), with the corresponding quiver velocities 
given by (see Eq. 20) 

<x>o 

fo j YOj (vo — £2 C j 
which can be rewritten as 



<2jB±.o 

mcko 


(46) 


v±oj =ajr]jce 


i(koz-coot ) 


(47) 
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with 


so that the dispersion relation is given by 


Oij = GjOL , 

0 

Tjj = -. 

fojYOj^O - n CJ 


(48) 

(49) 


The dispersion relation for the pump wave is (see Eq. 21) 


co 


2 

0 



COO 


fo j Yo j m Qc j 


(50) 


4.3 Parametric decays 


Now we introduce small perturbations on the zeroth-order 
solution given by the pump wave as described in Sect. 4.2. 
Thus, a physical quantity C is written in the form C = Co + 
8C , where Co corresponds to the pump wave solution, and 
8C « C. 

To first order in the perturbed quantities, we obtain the fol¬ 
lowing set of perturbed equations (Lopez et al., 2012) 


fo (iCi + Voj ' V ) t yoj8Vj + 8y J v oj) 

+ fo (Svj ■ V) (yojVoj) + Yoj + »o j • v j (»o. j&fj) 


— ( 8E + -vo/ x 8B ) + —8vj x Bo -— V8pj , 

m \ c J me mno 


d8 Pj 




-= me n n ■ -, 

3 1 0j 3 1 


3 8rij 


1 +« 0 V -Svj =0, 

3 1 

(51) 

V • 8E = 4n , 

(52) 


and 

4tt / x 1 3 8E 

v xSB = —J 2 < o( 8 n j v j +n j 8 v j) + ~~^f- (53) 

j 


We now introduce perturbations on the transverse quanti¬ 
ties (18) such that 8C± 8C Z <^C Z . We also assume 

that longitudinal and transverse perturbations have the form 


SC Z = Re ] 

= i [ce i(kz ~ M) + c*g- , '(***-‘"*0) 

. (54) 

and 



SC = c+ e« ( M-®+0 

+ c _ e i(k-Z-<o-t) ^ 

(55) 


respectively, where k+ = &o + k , k- = ko~k*, co+ = coo + co 
and co- = coo ~co*. With these definitions, and after a lengthy 
but straightforward calculation (see Lopez et al., 2012 for 
details), we obtain a set of equations for the variables r = 
(v+ e , v+ p , v^_ p , v Q , Up, b+, b*_ which can be put in 
matrix form as 


Ar = 0 , 


(56) 


F(k, co) = det(A) = 0. (57) 

The matrix elements A^ are given explicitly in Lopez et al. 
(2012). This turns out to yield a polynomial equation of de¬ 
gree 12 for co(k). 

In the cold, weakly relativistic case, we recover the results 
previously found by Munoz and Gomberoff (1998b). 

When there is no pump wave, a = 0, Eq. (57) yields the 
normal modes that propagate in the plasma. At this limit, it is 
not difficult to show that the dispersion relation (57) becomes 


F S d(Kcd) = D+D-SL = 0, 

(58) 

where 


S = —ft) 2 + v 2 k 2 , 

(59) 


2 CO 2 2/2 

L = —- — co + v:k , 

fo 


(60) 


and 


£>+ = coi 



( 61 ) 


D _ = 2 - c 2 k* 


fo 


( oO 

U-t 


^ \ 

w -+t/ 


and where we have defined 


(62) 


v s = c 


1 / , (mo)mo 
fo l + f'(m)Yk 


(63) 


Equation (59) suggests that v s can be defined as an effective 
sound velocity. 

In Lig. 11 we compare the effective sound velocity as 
given by Eq. (63) with the classical sound velocity, v s = 
sfkjfTfm. For low temperatures both results are equivalent, 
whereas for high temperatures the effective sound velocity 
tends to c/yf 3, which is a known result for hot relativistic 
plasmas (Landau and Lifshitz, 1959). 

Equation (58) shows that the normal modes correspond to 
an electroacoustic wave (given by S = 0), Langmuir waves 
(L = 0), and sideband waves ( D± =0). In the presence of 
the pump wave (a / 0) these normal modes couple and give 
rise to the parametric decays. We will study this in detail in 
Sect. 4.4. 


4.4 Numerical analysis of the dispersion relation 

We now solve the dispersion relation Eq. (57) numerically. 
To this end, for given values of a and /r, we choose a 
value of the normalized wave number of the pump wave, 
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Fig. 11. Normalized sound velocity as a function of normalized tem¬ 
perature =kftT/me 2 . Dashed line: classical result (/ = 1), 

v s /c = Solid line: relativistic result (see Eq. 63). The hori¬ 

zontal dotted line represents the limit value 1 / y/3. 

yo =koc/Q c , and we solve the dispersion relation of the 
pump wave to get the normalized pump wave frequency 
xq = coq/Q c , as outlined in Sect. 2.1 and Asenjo et al. (2009). 
Notice that, in particular, this implies a choice of a branch 
for the pump wave, since a given ko may correspond to up to 
three values of coq. 

Then we solve Eq. (57). In Fig. 12 the pump wave ( jco , 
yo) has been chosen to lie in the Alfven branch, near the ori¬ 
gin. Figure 12a shows the solutions of Eq. (57) for yo = 1, 
cop/Q c = 1 , 1 hi — 0.01 and a = 0 , which is equivalent to 
solving Eq. (58). There are 12 lines, two of which corre¬ 
spond to the electroacoustic modes that are the real solu¬ 
tions of S = 0. They are labeled as S+ and S- , corresponding 
to forward and backward propagating modes, respectively. 
Two other lines correspond to the Langmuir modes given 
by L = 0, and labeled as L + (forward) and L_ (backward). 
The other eight lines correspond to the eight real solutions 
of D± = 0. Four of these lines are parabolic and correspond 
to the electromagnetic branches in Fig. 1, and they are la¬ 
beled as D± in Fig. 12a. The two lines p± also correspond 
to solutions of D± , but they resonate at the sideband fre¬ 
quency co± = £ 2 C //o, which in the cold case is the positron 
gyrofrequency (see Eqs. 61 and 62). In terms of the normal¬ 
ized longitudinal frequency x = co/Q c , the resonance occurs 
at x = ±(— xo + l//o), as shown in Fig. 12. Similarly, e± 
are the two solutions that resonate at the sideband frequency 
co± = —Q c /fo (which, in the cold limit, is the electron gy¬ 
rofrequency), that is, at x = =b(—xo — l//o). 

We can see from Fig. 12a that there are several possible 
crossings between solutions of the dispersion relation. At 
these crossings, complex solutions can appear when a ^ 0 . 
Since the polynomial being solved has real coefficients, these 
solutions always occur as complex conjugate pairs, thus one 
of them has a positive imaginary frequency. Therefore, the 





Fig. 12. (Color online) Solution of the dispersion relation Eq. (57). 
Normalized wave number y = kc/£2 c vs. normalized frequency 
x = co/Q c for yo = 1, <Wp/ Q c = U l//x = 0.01. We take the pump 
wave at the Alfven branch, (a) a = 0. (b) Sector of interest for 
a = 0, (c) a = 0.1. Dotted lines represent the real part of the so¬ 
lution when it is complex. 

disappearance of real solutions when a ^ 0 implies the pres¬ 
ence of unstable waves, corresponding to the parametric de¬ 
cays of the pump wave. 
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Only the crossings involving modes satisfying energy con¬ 
servation co\ + C 02 = Icoo (Jt e Z + ) can give rise to wave cou¬ 
pling when the pump wave is switched on (a / 0). They are 
marked with the labels from (1) to (5) in Fig. 12a. Crossing 
(1) is actually not considered here because it is unstable only 
in a narrow range of temperatures (around /x -1 ~ 0.01) and 
its maximum growth rate is too small (10 -3 ) compared with 
the other crossings. Thus we focus our attention on cross¬ 
ings (2)-(5) that are displayed in detail in Fig. 12b for a = 0. 
Crossing (3) corresponds to a resonant coupling {t— 1), and 
crossings (2), (4) and (5) correspond to nonresonant cou¬ 
plings. 

Now we turn the pump wave on by considering a ^ 0. In 
Fig. 12c cy = 0.1, and we notice that some crossings become 
gaps. This means that at these crossings we have complex 
solutions whose real parts are indicated as dotted lines, while 
the real solutions correspond to the continuous lines. Hence, 
we now have instabilities, indicating wave coupling. For ex¬ 
ample, crossing (3) gives rise to a wave coupling between 
(S+, /?_), corresponding to an ordinary decay instability in 
which the pump wave decays into a forward-propagating 
electroacoustic mode of frequency co and a sideband wave 
of frequency co -. Crossing (4) leads to a nonresonant wave 
coupling between (/?+,/?_), in which the pump wave decays 
into two sideband waves of frequency co- and co +. In the case 
of crossing (5), there is no gap, and therefore no instability. 
Finally, crossing (2) represents an intersection of (/?+, e _) at 
the origin, so it corresponds to an electromagnetic modula- 
tional instability (Longtin and O. Sonnerup, 1986). 

As mentioned above, the appearance of gaps means that 
there exist modes whose frequency has a nonzero imaginary 
part. There are several possible decays simultaneously, but 
not all of them will be equally relevant. We can quantify this 
by calculating, for a given value of a , the maximum imag¬ 
inary part at each gap. By comparing these values for all 
gaps, we may then gain intuition regarding the relative im¬ 
portance of each decay for various values of a. This is shown 
in Fig. 13. We define the imaginary part of the frequency as 
T = Im(x), and in each gap we seek its maximum value r max 
(maximum growth rate), and we plot it against a and l//z. 

The first interesting feature we observe is that the resonant 
coupling (3) exhibits the expected linear dependence with 
a , for small a , while the nonresonant ones usually show a 
quadratic behavior for small values of a. Let us also note 
that, for every pump wave amplitude, the maximum growth 
rate always occurs for /i~ l < 0.3. In Fig. 13a, b and d we see 
that for /x _1 >0.1 the instabilities disappear, but Fig. 13c 
shows that crossing (4) is still unstable, so for high tempera¬ 
tures crossing (4) is the dominant one. 

Another interesting feature is that, for some crossings, 
there is a threshold value of a below which there is no in¬ 
stability, and that this threshold depends on the tempera¬ 
ture. This can be seen in Fig. 14a, which shows that cross¬ 
ing (3) leads to an instability only if a > 0.07, approximately, 
when l//x = 0.03, whereas there is no threshold for the larger 



1 /p 


Fig. 13. (Color online) Maximum growth rate for various pump 
wave amplitudes a and temperatures l//x. Notice that the scale 
is logarithmic scale for temperature, (a) Crossing (2), (/?+, e_). 
(b) Crossing (3), (S+, p_). (c) Crossing (4), (p+, /?_). (d) Crossing 
(5),(*_, p_). 
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Fig. 14. Maximum growth rate vs. pump wave amplitude a. Nu¬ 
merical labels for crossings are taken from Fig. 12. (a) Crossing 
(3), (£+, p -). (b) Crossing (5), ( e _, p-). 




Fig. 15. Imaginary part of frequency vs. normalized wave number, 
for various values of the pump wave amplitude. Numerical labels 
for crossings are taken from Fig. 12. (a) 1 //x = 0.01. (b) 1 //x = 
0 . 02 . 


temperatures shown in the plot. In the case of crossing (5) 
(Fig. 14b), there is a threshold for all temperatures. When 
temperature is increased, smaller values of a, are needed to 
trigger the instability. 

In Fig. 15 we have plotted the imaginary part of the fre¬ 
quency as a function of the wave number, for different values 
of a and a fixed value of /x. This figure gives us an indication 
of the wave band over which the wave coupling occurs. We 
have labeled the instabilities as in Fig. 12a. In Fig. 15a we 
take l//x = 0.01, and as we expect for this range of values of 
a , the resonant coupling (3) is the dominant one, being larger 
and broader than the other, nonresonant, couplings (2), (4) 
and (5). This also occurs in Fig. 15b, for l//x = 0.02. How¬ 
ever, the wave band of active modes of crossing (3) has been 
significantly reduced, while the ones for crossings (2) and (4) 
have slightly increased. In Fig. 15, crossing (5) does not ap¬ 
pear because it is stable for these temperatures. 

Figure 15 suggests that the dominant coupling may be 
different for different values of a and /x. We can appreci¬ 
ate this in Fig. 16, where we plot the maximum growth rate 
among all crossings for various temperatures and pump wave 



1 /p 


Fig. 16. (Color online) Maximum growth rate for various pump 
wave amplitudes a and temperatures l//x, in logarithmic scale for 
temperature, considering all existent crossings. 


amplitudes. Discontinuities in the derivative of the plotted 
surface are an indication of a change of regime, where the 
maximum growth rate changes from one crossing to another. 
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Fig. 17. Maximum growth rate vs. pump wave wave number, for 
a = 0.04 and various temperatures. The vertical lines indicate the 
value of jo at which dco/dk = 0. 


In particular, for low temperatures the dominant behavior 
is given by crossing (3), while for higher temperatures it is 
given by crossing (4), which turns out to be the only cou¬ 
pling present. 

A similar analysis to the one presented in this section can 
be done for a pump wave in the electromagnetic branch and 
in the anomalous zone of the Alfven branch (see Lopez et al., 
2012 for full details), with qualitatively similar results. 

One interesting finding is that if the pump wave is in the 
Alfven branch, r max is about an order of magnitude larger 
than in Fig. 16, when the pump is a wave in the normal 
zone of the Alfven branch. This difference can be observed in 
Fig. 17, where we set a pump wave amplitude a = 0.04 and 
a temperature \i ~ x , and the maximum growth rate is plotted 
for various values of yo in the Alfven branch. When yo ap¬ 
proaches the anomalous zone, the instability becomes larger. 


5 Numerical simulations 

In this section we present preliminary results from relati¬ 
vistic particle simulations that are consistent with the find¬ 
ings in Sects. 2, 3, and 4. Simulations have been done with 
a relativistic PIC code, with periodic boundary conditions, 
and initial particle thermal velocities have been taken from a 
Maxwell-Boltzmann-Juttner distribution function. The de¬ 
tailed description of the code, as well as more complete re¬ 
sults, will be presented elsewhere (Lopez et al., 2014). 

First, by choosing an initial condition where particles only 
have thermal velocities, we can obtain the normal mode dis¬ 
persion relation from the power spectrum of magnetic fluc¬ 
tuations. This is shown in Fig. 18. 

A good agreement with both fluid (Sect. 2) and kinetic 
(Sect. 3) results is found, as shown in Fig. 19, where the 
curves corresponding to the appropriate dispersion relations 



Fig. 18. Power spectrum of transverse magnetic fluctuations from a 
PIC simulation, for a = 0.01, <x> p /£2 c = 1. (a) /x = 100, (b) /x = 7. 




y 


Fig. 19. Same as Fig. 18, but with curves from the analytical models 
superimposed. Blue dotted line: fluid dispersion relation, Eq. (21); 
black dotted line: kinetic dispersion relation, Eq. (40). (a) /x = 100, 
(b) ii = l. 
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Fig. 20. Power spectrum of density fluctuations for the PIC simula¬ 
tion. The pump wave satisfies Eq. (21), with a = 0.2 and /z = 100, 
and is taken from the Alfven branch, with yq = 0.49, xq = 0.27. 
(a) Simulation results; (b) lines correspond to the fluid dispersion 
relation of the parametric decays, Eq. (57), for the same parame¬ 
ters. 


(Figs. 1 and 6) have been superimposed. In particular, both 
expected branches are observed, the electromagnetic and the 
Alfven one. Regarding the Alfven branch, in the simulation 
results it cannot be distinguished for wave numbers beyond 
a certain critical value, which suggests that it is damped due 
to kinetic effects. 

Figure 19 also shows that the effective plasma frequency 
decreases with temperature (see Figs. 2 and 7), and that 
the Alfven velocity decreases with temperature as well (see 
Eq. 23 and Fig. 8). When the effective plasma frequency and 
Alfven velocity are measured from the simulation results, 
they turn out to be consistent with the fluid and kinetic mo¬ 
dels in Sects. 2 and 3. 

We have also studied the parametric decays of a circularly 
polarized pump wave, following the model in Sect. 4. Fig¬ 
ure 20a shows the power spectrum of density fluctuations, 
when the initial condition is a pump wave consistent with 
Eqs. (10) and (11), with a = 0.2 and /z = 100. Various modes 
corresponding to the parametric decays of the pump wave 
can be identified, with major wave activity near y ~ 0.8. In 
Fig. 20b the dispersion relation for the parametric decays in 
the fluid model for the same parameters has been superim¬ 




Fig. 21. Growth rate for the decay instability shown in Fig. 20. 
(a) Fluid result for the (S+, p-) crossing, (b) Simulation results 
for the same parameters. 

posed, showing the consistency of the simulation and ana¬ 
lytical results. In particular, the spectrum peak near y ~ 0.8 
can be identified as the result of the nonlinear coupling be¬ 
tween a sound wave 5+ and a p- sideband wave (see, e.g., 
crossing (3) in Fig. 12b). 

Finally, growth rates for the decays can be estimated from 
the simulation. This is shown in Fig. 21, where the fluid and 
simulation results are shown. Again, both models are consis¬ 
tent. 

6 Conclusions 

We have studied the propagation of large-amplitude, circu¬ 
larly polarized electromagnetic waves in electron-positron 
plasmas, using fluid and kinetic models, and PIC simulations. 
Relativistic effects in the particle velocities have been intro¬ 
duced both for the quiver velocity induced by the wave, and 
for the random velocities induced by temperature. In the case 
of the fluid model, relativistic temperatures have been consis¬ 
tently introduced through a factor / related to the enthalpy 
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density. For the kinetic model and the simulations, relativistic 
temperatures have been introduced by assuming that the par¬ 
ticles follow a Maxwell-Juttner distribution function. Prop¬ 
agation of the wave has been assumed to occur along a con¬ 
stant background magnetic field. 

In the case of the fluid model, an exact solution of the re¬ 
lativistic fluid equations can be found. A critical assumption 
is the fact that particle velocities are purely transverse with 
respect to the background magnetic field, which is true for 
circularly polarized electromagnetic waves of arbitrary am¬ 
plitude. In particular, this means that no pressure or density 
fluctuations appear, and that / is constant. We plan to explore 
the possibility of extending the above analysis to include a 
pressure tensor (Shadwick et al., 2004), a case in which we 
do not expect to find exact solutions, but which can certainly 
be studied numerically. 

The fluid results show that two branches are present, an 
electromagnetic and an Alfven branch. The electromagnetic 
branch has a lower cutoff at an effective plasma frequency 
cof, which decreases with temperature and wave amplitude 
due to the relativistic increase in effective mass of the parti¬ 
cles. 

The Alfven branch, on the other hand, shows several in¬ 
teresting features. For very small wave numbers it is non- 
dispersive, with an Alfven velocity va that depends only on 
temperature, not wave amplitude. As wave number increases, 
a maximum frequency co cr i t , corresponding to a critical wave 
number /c cr j t , is reached. This upper cutoff depends both on 
wave amplitude and temperature. Thus, as temperature in¬ 
creases, Alfven waves are confined to an ever narrower fre¬ 
quency bandwidth. 

In a fluid model, important effects such as Landau damp¬ 
ing are absent. Thus, a kinetic description is necessary. We 
have considered this problem too, studying the same prob¬ 
lem, but based on the relativistic Vlasov equation, in order 
to account for kinetic effects in the direction of the magnetic 
field. The dispersion relation of the waves is found and nu¬ 
merically solved. Results are qualitatively consistent with the 
fluid model: the effective plasma frequency decreases, and 
the frequency range of the Alfven branch is reduced when 
the temperature is increased, so that electromagnetic waves 
propagate in a larger frequency range with respect to the cold 
case, whereas for Alfven waves the opposite occurs, being 
confined to a frequency range that is narrower than in the 
cold case. 

However, quantitative results are different in both models. 
Regarding the electromagnetic branch, it is found that the ef¬ 
fective plasma frequency is always larger in the kinetic model 
as compared with the fluid one. Thus, kinetic effects reduce 
the transparency of the plasma. As to the Alfven branch, in 
both models it is found that the Alfven speed decreases with 
temperature. In the fluid model, the Alfven speed is zero for 
infinite temperature. However, in the kinetic case there is a 
critical, nonzero value of the temperature at which the Alfven 
speed is zero. For temperatures above this critical value, the 


Alfven branch is suppressed. The critical value depends on 
the background magnetic field. For larger values of the mag¬ 
netic field, the critical value increases, and thus Alfven waves 
can propagate for larger temperatures. 

We also studied the fluid dispersion relation of the waves 
for nonzero drift, thus generalizing the numerical results in 
Sect. 2. Besides, this analysis allowed us to study the depen¬ 
dence of the relativistic Lorentz factors on the drift, which 
was later used to study the kinetic model as detailed in 
Sect. 3.4. 

We are also interested in the nonlinear evolution of the 
large-amplitude waves studied in Sects. 2 and 3. Thus, we 
have also studied the parametric decays of this wave, in the 
context of the fluid model. We study the appearance, the 
growth rate and instability range of the decays for various 
temperatures and pump wave amplitudes. We also consider 
pump waves in three possible regions: in the Alfven branch 
(near the origin and in the anomalous zone), and in the elec¬ 
tromagnetic branch. There are both resonant and nonreso¬ 
nant couplings when the pump wave is present. For high 
enough temperatures, though, most of the wave couplings 
disappear. For certain wave couplings there is a threshold 
pump wave amplitude, such that instability arises only above 
the threshold. 

Finally, we have shown the consistency of the results 
above with those obtained from a relativistic PIC simulation. 
The existence of an electromagnetic and an Alfven branch, 
the decrease in the effective plasma frequency and the Alfven 
velocity with temperature, the confinement of the Alfven 
branch to a narrow region of wave numbers and frequencies 
as temperature is increased, and the dispersion relation of the 
parametric decays, have been shown to agree qualitatively 
and in some cases quantitatively with both analytical models. 

We expect to improve further on these results. For in¬ 
stance, we have not studied the damping due to kinetic 
effects. Another improvement would be to consider a ki¬ 
netic model with a three-dimensional, rather than one¬ 
dimensional, particle velocity distribution. Also, kinetic ef¬ 
fects on the parametric decays need to be introduced, and 
some analytical results are available for the relativistic case 
(Munoz et al., 2006). We have also studied the nonlinear evo¬ 
lution of electromagnetic waves in this system by means of 
the NLS equation, finding conditions for the existence of en¬ 
velope solitons (Asenjo et al., 2012), localized structures that 
can also be found in our simulation runs, and we plan to in¬ 
vestigate further these and other issues related to the prop¬ 
agation of nonlinear waves in relativistic pair plasmas else¬ 
where. 
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